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1. Introduction 

A quantitative understanding of interesting properties of the flavor-singlet mesons is an impor- 
tant subject in lattice QCD. The famous U{1) problem is a long-standing issue albeit past ceaseless 
efforts. It is also well-known that the flavor-singlet and octet vector mesons mix with each other 
almost ideally, though this has not yet been confirmed from first principles. 

In this article, we report on our study of the flavor-singlet pseudoscalar (PS) and vector mesons. 
There are two salient features of this work: i) we simulate Nf = 2+1 QCD including the effects of 
dynamical strange quarks, which have been often ignored in previous studies of the flavor-singlet 
mesons, and ii) we use the all-to-all quark propagator [|l]] to calculate disconnected meson correla- 
tors, which induce the meson mixings and the mass splittings from flavor-non-singlet mesons. 



2. Simulation method 

Our gauge configurations of Nf = 2+1 QCD are generated on a 16^ x 48 lattice using the 
Iwasaki gauge action and the overlap quark action. We also introduce a topology fixing term [§] 
into our lattice action to reduce the computational cost, and simulate only the trivial topological 
sector 2 = at this stage. The lattice spacing determined from Fj^ is 0.100(5) fm. We take four 
values of the degenerate up and down quark masses mi = 0.015, 0.025, 0.035 and 0.050, which 
cover a range of the pion mass from 350 to 610 MeV. Two values = 0.080 and 0.100 are chosen 
for the strange quark mass. The physical quark masses fixed from Mji and Mk are m/ phys = 0.002 
and phys = 0.065. Statistics are 2,500 HMC trajectories at each combination of nii and m^.. We 
refer readers to Ref.[Q] for further details on our gauge configurations. 

We measure PS and vector meson correlators using the all-to-all quark propagator. For each 
configuration, we prepare 160 low-lying modes \ u^^^) {k=l,...,Ne{= 160)) of the overlap- 
Dirac operator D{m), where m is the valence quark mass. Their contribution to the quark propagator 
is calculated exactly. The higher modes are taken into account stochastically by the noise method. 
We prepare a single noise vector rj for each configuration, and dilute [|l|] it into Nd = 3 x 4 x Nt/2 
vectors ri^'^^ (d = \,..,Nd), which have nonzero elements for a single combination of color and 
spinor indices and at two consecutive time-slices. The all-to-all propagator can be expressed as 

Dim)-' = £ vi? w«t (^^, = (2. 1) 
k=l 

with two set of vectors 

v« = (^,...,^,x2\...,xif4, w« = {«(^...,«W,r7(»,...,r7(-^)}, (2.2) 

where Xm'^ is the solution of the linear equation 

D(m)xif = (1-£m«m«^)t]W. (2.3) 

k 

We then construct the following meson field at the temporal coordinate t with the Dirac matrix 
r and a smearing function (r) 

('«;0 = E0(r)w(x + r,O('^^rvm(x,O('). (2.4) 
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Figure 1: Connected (left diagram) and disconnected meson correlators (right diagram). We denote the 
flavor matrix and the smearing function for the sink meson operator by Fa and 0,, and those for the source 
operator by F/, and (j)j. The Dirac matrix is r= 75 for PS mesons, whereas we average over r= 71 2.3 for 
vector mesons. The masses of propagating quarks are denoted by m^(,) . 

The connected and disconnected correlators shown in Fig. |l| can be calculated from these meson 
fields. The PS meson correlators, for instance, are given by 

Cp,ab,ijiAt) = I I (F,),,,(f,),,,,^J:'^)K;f + A0^J;2.('«,';0, (2-5) 

' t q'=u,d,sl=\ q=ii,d,sk=\ 

For simplicity, we often suppress the indices of the smearing functions (/ and j) in the following. 

In this study, we consider the PS and vector mesons in two different flavor bases: i) light and 
strange mesons, P/ ^ and V/,s , with their flavor matrices F/ = (l/\/2)diag[l, 1,0] andFi.=diag[0,0, 1], 
and ii) octet and singlet mesons, P% q and o, with U (3) generators o for the flavor matrices {F% = 
and Fq = Tq). We refer to these bases as the light-strange and octet-singlet bases, respectively. 
The full correlator of these mesons including the disconnected contribution is given by 

G{py}^ab,ij{M) = C{py},ab,ij{At) - D{py}^ab,iji^t) {a,b G {l,s} ora,b G {8,0}). (2.7) 

We calculate all possible correlators with the following five different choices of the smearing func- 
tion (namely /, j = 0, ...,4) 

</>o(r) = 5r,o, 0i(r) °< exp[-0.4|r|], ^2(r) °< |r| exp[-0.4|r|], 

(z.o) 

(^3(r) oc exp[— 1.0|r|], 04(r) = constant 

with the normalization Lr I^K*")!^ = 1- The calculation of all these meson correlators is computa- 
tionally cheap, once we prepare the v and w vectors of Eq. (^^. 

Since the quark propagator is decomposed into low- and high-mode contributions, the discon- 
nected correlators can be divided into four contributions, i.e. D = D^^ + D^^ + D^^ + Z)™. We 
calculate these four contributions separately in our measurement. 

3. Meson correlators 

In Fig. we show an example of the light PS and vector meson correlators. We observe that, 
at relatively small At, the disconnected piece Dpu in the full correlator Gpji is not large and it is 
dominated by the low-mode contribution £)p^/. Therefore, we may safely ignore the high mode 
contributions to Dp//, namely D^i, Dpfj and Dp^, to calculate the full correlator as 

G?;l-jj{^t) = CpahJjiAt) - D|;^,/,,/y(AO. (3.1) 
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Figure 2: Light PS (left panel) and vector meson correlators (right panel) at (m/,ms) = (0.025,0.080) with 
exponential smearing 0i for source and sink. Open circles, squares and triangles show the full (G{/)v},//)> 
connected (C^pyjji) and disconnected correlators (D^pyj n), respectively. We also plot the low-mode con- 
tribution to the disconnected piece D^pyj n by diamonds, and the full correlator G^jTy^n defined in Eq .^.l[ ) 
by filled circles. We note that the full PS meson correlator has a constant term as discussed in Section ^. 



Figure ^ actually shows that Gp,ii is well approximated by Gpj^^^ in the whole region of At. 

As shown in the same figure, the vector meson full correlator Gyji turns out to be noisy at 
relatively large At. The large uncertainty mainly comes from those of high-mode contributions 
Dyii ' ' due to the noise method with the small number of noise samples. We observe that 
Gv^ii at small At is well approximated by Gyjj^^ defined as in Eq. (^?T|), and expect that the high- 

f HL HH \ 

mode contributions Dyu ' remain to be small at larger At since they mainly describe short 
distance physics. Then Gyji is expected to be well approximated by GyJ;^^ also at large At. 

From these observations, we use the meson correlators G^jTy^^f^ ignoring the noisy contribu- 
tions D^^p^'yl^if^^ to study the spectrum of the fiavor-singlet mesons. The superscript ''D=LV' is 
suppressed in the following for simplicity. 



4. Vector mesons 

We plot the vector meson correlators in Fig. ||. In the light-strange basis, the off-diagonal 
correlators Gy^^is^si] are about two orders of magnitude smaller than the diagonal ones Gyj^u ss}- 
There is no such large hierarchy in Gy {gg oo,80,08} i^i the octet-singlet basis. Since the off-diagonal 
correlators induce the meson mixing, the above observations indicate that the mixing of the vector 
mesons is close to the ideal mixing: namely, V% and Vq mesons mix significantly with each other to 
form ft) and mesons, which are well approximated by V/ and Vg. 

For a more quantitative examination, we solve the generalized eigenvalue problem (GEVP) 

C{At')-^l'^C{At)C{At')-^l'^Un = XnUn (?i = 0, 1), (4.1) 
where C(A?) is 2 x 2 correlator matrix with specified smearing functions for source and sink 

^ ' \ Gy,08(A0 Gy,oo(AO 
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Figure 3: Vector meson correlators in the light-strange (left panel) and octet-singlet bases (right panel). 
Both panels show results with the exponential smearing function at (m;,ms) = (0.025,0.080). Filled and 
open symbols represent diagonal and off-diagonal correlators, respectively. 



The creation operators of the energy eigenstates, namely, <p and co mesons, are determined from 
the eigenvectors m„. We obtain the following relation for the local operators 



= 0.84(5)^8- 0.55(7) Vb 
(0 = 0.55(7)^8+ 0.84(5) Vo 



1.00(1) V,- 0.04(9) V/ 
0.04(9)V, + 1.00(1)V, 



(4.3) 



which implies the ideal mixing of the vector mesons. We observe that this relation of vector meson 
operators does not change significantly with other choices of the smearing functions. 

In Fig. ^ we plot the effective masses of co and ^ mesons determined from the eigenvalues A„ = 
exp[—E„{At — At')]. By taking a sufficiently large At', the effective masses show small dependence 
on At as well as on the smearing functions. The same figure also shows that Mq, and are close to 
those from the connected correlators of the light and strange mesons Cv^{iijs}- This is because the 
diagonal disconnected correlators {// ^^.j are small as seen in Fig. ^ and hence they have small 
effects to Ma, and M^. 

In this analysis, we extrapolate M^(0) to the physical point using a simple linear form 
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Figure 4: Left panel: effective masses of o (solid symbols) and mesons (shaded symbols) at (w^rf, Wj) = 
(0.025,0.080). We also plot effective masses of Cyu (open circles) and Cv,^ (open squares). Right panel: 
chiral extrapolation of . Dotted and dashed lines show fit lines at m^ and phys, respectively. 
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As shown in Fig. ^ this fit describes our data reasonably well with ;f^/d.o.f ~ 1.3. We obtain 
Mco = 909(25)stat(-98)sys MeV andM^ = 1102(13)stat(-97)sys MeV, where the systematic error 
is estimated by including a higher order term <io(0)'n;^''^ |Q] and by using a different input Mq, to fix 
the lattice spacing. These results are consistent with the experimental values Mg, = 783 MeV and 
= 1019 MeV. Note, however, that our data may suffer from significant finite volume corrections 
at two smallest quark masses m/ = 0.015 and 0.025, where 2.8 <M;r^ < 3.2. We are planning to 
extend this work to a larger volume 24^ x 48 for a more precise comparison with the experiment. 



5. PS mesons 



Figure |5| shows correlators of the light and strange PS mesons. In contrast to the vector mesons, 
the off-diagonal correlators Gpg^ sij are not so small compared to the diagonal ones Gp{// „| in the 
light-strange basis. This leads to a significant strange (light) quark component in rj (rj'). For the 
local operators, we obtain 



T] =0.96(l)P/-0.28(3)/'v, T]' = 0.28(3)P/-0.96(1)P,. 



(5.1) 



As in phenomenological analyses more unambiguous determination of the mixing matrix could 
be provided by constructing the local Tj and rj' operators so that their decay constants reproduce 
the experimental values. We leave this for a future study. 

As predicted analytically ^ and as seen 
in Figs. ^ and disconnected contribution Dp^ah 
induces a constant term in the full correlator 
Gp^ab at fixed topology, e.g. 



miGpji 



Ar- 



V 



1 



(5.2) 




which is suppressed by l/V. While this term 
is useful to determine the topological suscepti- 
bility Xt [Q]> this forces us to use Gpat at small 
At to extract the PS meson masses M,, and M'^ . Figure 5: Light and strange PS meson coiTelators. 

To eliminate excited state contamination at such small At, we solve the GEVP with the 10 x 10 
correlator matrix including the smearing degrees of freedom 



/ GpMfloi^t) ■■■ Gp3o,04(A?) 



C{At) 



(5.3) 



\Gpfl^^4o{At) ••• Gpfl0A4{At) 
Effective masses of rj and rj' mesons are plotted in Fig. |^ Although rj seems to be lighter than rj', 



the existence of the constant term in Eq. ( |5.2| ) leads to a large uncertainty of M^/: 10-15% already 
at simulated quark masses. 

From a linear chiral extrapolation in terms of mi and nis, we obtain Mjj =639(50)stat MeV and 
Mjj/ = 840(136)stat MeV at the physical point. These are consistent with the experimental values 
Mrj =548 MeV and M^/ = 958 MeV, though the statistical significance of the rj' -rj mass spUtting 
(1.4 a) is not sufficient. 
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6. Conclusion 
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In this article, we report on our study of 
the flavor-singlet mesons in Nf = 2 + 1 lattice 
QCD using the all-to-all quark propagator to 
calculate the disconnected correlators. For 
the vector meson, we observe that the small 
disconnected contributions in the light-strange 
basis lead to the almost ideal mixing and the 
small mass shift. This is consistent with the 
experimental fact Ma, ~Mp, and explains why 
the previous calculations of ignoring the 
disconnected contributions show reasonable 
agreement with experiment [^. 

We need to improve the accuracy of M^i 
to establish the rj'-rj mass splitting. This 
could be done by simulating non-trivial topo- 
logical sectors as well as by suppressing the fixed topology effects in Eq. (5^) on a larger lattice. 
The latter is also important to suppress finite volume corrections to the PS and vector meson masses 
for a more detailed comparison with the experiment. Such simulations on a 24^ x 48 lattice are in 
progress. 



Figure 6: Effective masses of rj' (top panel) and rj 
(bottom panel) at (m„rf,m,s) = (0.025, 0.080) with differ- 
ent choices of At'. For a comparison, we also plot effec- 
tive mass from the connected strange correlator Cp^ss- 
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